Longitudinal Force on a Moving Potential 
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We show a formal result of the longitudinal force acting on a moving potential. The potential 
can be velocity-dependent, which appears in various interesting physical systems, such as electrons 
in the presence of a magnetic flux-line, or phonons scattering off a moving vortex. By using the 
phase-shift analysis, we are able to show the equivalence between the adiabatic perturbation theory 
and the kinetic theory for the longitudinal force in the dilute gas limit. 

PACS numbers: 67.40.Vs 



I. INTRODUCTION 



In recent work on the theory of the tpajisverse force on 
vortex moving relative to a superfluidtJcl we have used 
a technique in which the force is calculated by apply- 
ing a pinning potential which is made to move, and the 
reaction on the pinning potential is calculated pertur- 
batively. A number of people have objected privately 
that this method does not give the longitudinal force 
correctly. We do not agree with this objection, since 
we think that it can be used to determine the longitudi- 
nal force. For simple problems our method gives results 
that are in agreement with those obtained by more fa- 
miliar methods. This paper examines this question for 
some special cases. We conclude that any differences 
from other perturbative methods which may occur are 
due to uncertainties in limiting procedures. 

For the component of force transverse to the vortex 
velocity on a vortex in a superfluid Thouless, Ao and 
NiuEl (which we refer to as TAN) found that its value 
is independent of the pinning potential with which the 
vortex is made to move. The magnitude of this trans- 
verse force per unit length is equal to the circulation 
of momentum density at large distances multiplied by 
the velocity of the vortex relative to the superfluid, to 
lowest order in the vortex velocity. In ideal supercon- 
ductors a similar method gave a transverse force equal 
to the product of the vortex velocity, and the circula- 
tion of canonical momentum densitya. However, this 
result is controversialjlj. The discrepancy may come 
from using different calculation schemes, or it may come 
from inadequate handling oi boundary conditions in some 
or all of the calculationsH. TAN's derivation uses the 
time-dependent perturbation theory in contrast to the 
usual approach using scattering theory and kinetic equa- 
tions. It is known that the kinetic method and the 
linear-response theory do not have the same domains of 
validity^, which motivates us to examine the applicability 
of the force formula more carefully. 

In Sec. II we consider the effect of a moving short- 
ranged spherically symmetric potential on a system of 
nonintcracting particles, and show that the longitudinal 



force obtained from the method of TAN is, to lowest or- 
der in the velocity, identical to that which is obtained 
from kinetic theory. In Sec. Ill we consider the long- 
ranged Aharonov-Bohm type potential, and again find 
no difference in the longitudinal force. The transverse 
force seems to be more sensitive to limiting procedures. 



II. THE LONGITUDINAL FORCE 

The Hamiltonian describing a system of noninteract- 
ing particles in a system which is uniform apart from a 
spherically symmetric scattering potential whose center 
is time-dependent has the form 
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It is convenient to work with a velocity which is switched 
on adiabatically, so that ro = vy e 7 * . Since we will only 
be considering effects linear in vy the frequency depen- 
dence of the response can be obtained by analytical con- 
tinuation of the 7 dependence into the complex plane. 
The method used by TAN involves calculating the expec- 
tation value of the force on the potential to first order in 
vy, using the instantaneous eigenstates of H as a basis. 
This gives 
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vy • V + h.c. 



where l^a) is an instantaneous eigenvalue of H with 
eigenvalue E a , f a is the equilibrium probability that the 
state is occupied, and V a is the projection operator off 
the state a. The first term is not proportional to the vor- 
tex velocity Vy, and vanishes if the system has spherical 
symmetry. We can simplify Eq. (0) further by substitut- 
ing Vol^o) with V H, which is V V, because V H is 
the commutator of Vo and H , 



V |* a ) = 
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This gives the resultant longitudinal force as 
fa - fp 



E n — En 
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In the limit that 7 goes to zero, which is the case that 
we consider in the rest of this paper, the last term gives 
a delta function of the energy, so that the first factor on 
the right hand side gives df/dE. 

For noninteracting particles the sum over states can 
be replaced by a sum over single-particle wave functions, 
with the f a replaced by fermionic or bosonic occupa- 
tion probabilities. For a central potential the energy- 
conserving matrix elements in Eq. (Q) can be expressed 
in terms of phase shifts, using the relation between the 
plane wave and spherical wave expansions, 
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where V is the system volume, and the radial wave func- 
tions at large distances are 



Ai{kr) = e lS ' [cos Si ji(kr) — sin 6ini(kr)] 



(6) 



The matrix element which needs to be evaluated for the 
longitudinal force given by Eq. (^) is 
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if Tq is taken to be in the z direction. The last angular 
integral can be carried out by using the following recur- 
rence relation, 



cos 6 Y, 
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where the orbital angular momenta are changed by one 
because the gradient is a tensor operator of rank one. 
Taking the square of its modulus and integrating k a and 
kfj over all possible directions will get rid of the remain- 
ing spherical harmonics, and end up with a summation 
over orbital angular momenta, 
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The summation over magnetic quantum numbers can be 
easily carried out because there is only one independent 
direction left after the integration. The radial integral 
gives the difference between phase shifts (see Appendix 
A), 



2m J drr 2 A^ +1 (kr)^Ai(k7 
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There is no contribution to the force from bound states 
because their wave functions vanish at large distances, 
and resonance states may potentially play an important 
role here. So the longitudinal force can be written as an 
integral over all scattering states, 
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This result was obtained in a similar way by Bonig and 
Schonhammeiu. 

There is a very similar calculation for scattering by an 
axially symmetric potential, for which the wave functions 
are 



im(tj>-tj>- k ) 



(12) 



The calculation for the matrix element is almost the same 
except for a different normalization constant, 
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After integrating all incident directions, only the radial 
integral remains, 
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and this gives the connection to the phase shifts in the 
same fashion, 

[ , .* dV A 
2m J dppA m+1 —A m 
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This gives the longitudinal force as 
p r at. u2 
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where L is the length of the system in the axial direction. 

Equations ([ll]) and ( |l6|) are identical to the equations 
obtained from the standard kinetic theory argument. In 
kinetic theory the usual procedure is to keep the scat- 
tering potential fixed and to shift the equilibrium distri- 
bution so that particles have an average momentum at 
large distances given by an undetermined multiplier vy . 
The force is the momentum transfer per unit time from 
particles to the scattering center. In d dimensions this is 
related to the differential scattering cross section er(k, 0) 
by 



F = Vy 



d d k 



(2tt) 



\ N f It 

m v v ) / k k- 

m 



m Vy 



dO(l 



m Vy, I 



Nv v 
dVL 



dkk 



d+i 9fk 
dk 



(17) 



(18) 



k 2 
COS O k 



6 k / dfi(l- cos 6)a{k,6) , 



(2?r) d 

where, on the right side of the equation, only the first 
order term in the velocity has been kept. The angular 
integral over all incident directions can be carried out 
independently because of the isotropy, and gives 

dtt k cos 2 9 k = ^tttt-t-, . (19) 



T(d/2)d 

The last scattering cross section integral is often called 
the transport cross section, a tr . For a spherically sym- 
metric potential in three dimensions the scattering am- 
plitude can be expanded in Legendre polynomials as 

1 
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where 8i is the phase shift for the l-th partial wave. The 
differential cross section is just the modulus squared of 
the scattering amplitude, and the transport cross section 
is 
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This is in agreement with Eq. (|1 lj) . 

In two dimensions, the scattering amplitude is ex- 
panded in Fourier series, 

n<t>) = &i5m sin<5 ™ ^ > ( 22 ) 

m 

and then transport cross section follows similarly, 
atr = / d(j}(l — cos 4>)a(4>) 



= r E sin2 (<Wi - S m ) , 

m 

in agreement with Eq. (Jl6[). 



(23) 



III. LONG-RANGED POTENTIALS 



For problems guoh as the scattering of phonons by a 
quantized vortexBtL9 or the scattering of charged particles 
by an Aharonov-Bohm flux lineElrEj we need to consider 
potentials which fall off with distance like 1/r 2 and for 
which the phase shifts have a nonzero limit for large an- 
gular momentum. A potential of the form 



2m V(p, fy) = - 2ia-^ + fi 



(24) 



This gives an approximate form for scattering of phonons 
from a quantized vortex for a = ltKo/2irc, and for scat- 
tering of charged particles from an Aharonov-Bohm flux 



line for (3 
function 



where 



of. The radial wave function is the Bessel 
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v(m) = y to 2 + 2am + (3 , 
so that the phase shifts have the form 
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In place of Eq. (|l4|) we have 
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The integral on the right side of this equation can be eval- 
uated in terms of standard integrals of Bessel functionsEJ, 
and is equal to 



(k/n) cos{[i/(m + l) - v{m)]ir/2}e i ( s ™- Sm+l) 

2a[v(m + l) 2 - v(m) 2 - 1] - 4(am + (3) 
[v{m + l) 2 - v(m) 2 ] 2 - 2v{m + l) 2 - 2v(m) 2 + 1 



(29) 



Substitution of Eq. (|26|) into the fraction shows that it 
is equal to unity, while Eq. (^7|) shows that the cosine 
is equal to sin(5 m+ i — S m ) so that the force is given in 
exactly the same form as in Eqs. ( [Tif ) and (|2^). 

For this type of potential the forward scattering am- 
plitude diverges like 1/0, but for the transport cross sec- 
tion the factor of 1 — cos <f> cancels the divergence. For 
the transverse cross section the situation is much more 
delicate, as is-ishown in the recent discussion by Wexler 
and ThoulessEl. 
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IV. SUMMARY 

We have shown that, both for a short-ranged central 
potential, and for a long-ranged potential of the sort 
given by a flux line or vortex, calculation of the longi- 
tudinal force on a moving potential in a noninteracting 
background by kinetic theory and by time-dependent per- 
turbation theory give the same results to lowest order in 
the velocity. 

We thank Ping Ao for helpful correspondence. This 
work was partially support by NSF Grant No. DMR- 
9528345. 



APPENDIX A: RADIAL SCHRODINGER 
EQUATION 

Here is the key identity for establishing the equiva- 
lence between the two theories. Since we did not specify 
a particular pinning potential, any general property of 
the matrix element could be derived from the radial wave 
equation which in three dimensions is 

l^[rMkr)] + \k*-2 mo V(r)- 1 -^ 



Mkr)=Q. 
(Al) 



A general relation between different Us can be derived as 
following, 
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To evaluate the radial part of the matrix element in Eq. 
(0), we first integrate it by parts, replace V(r)AjJkr) us- 
ing the radial wave equation, and then use Eq. ( |A2| ) with 
n = 1, 
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Using the asymptotic expansion for the Bessel func- 
tions, the boundary term gives the difference between 
two phase shifts. Similar relation could be established 
from the wave equation in two dimensions, or could be 



done by directly substituting AAx) — > w — A m (x) with 
l = m- 1/2, 
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